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Abstract 

After giving an overview of the history and idea of neutron stars, I shall 
discuss, in part one of my lectures, the physics of dense neutron star matter. 
In this connection concepts like chemical equilibrium and electric charge neu- 
trality will be outlined, and the baryonic and mesonic degrees of freedom of 
neutron star matter as well as the transition of confined hadronic matter into 
quark matter at supernuclear densities are studied. Special emphasis is put 
onto the possible absolute stability of strange quark matter. Finally models for 
the equation of state of dense neutron star matter, which account for various 
physical possibilities concerning the unknown behavior of superdense neutron 
star matter, will be derived in the framework of relativistic nuclear field the- 
ory. Part two of my lectures deals with the construction of models of static 
as well as rapidly rotating neutron stars and an investigation of the cooling 
behavior of such objects. Both is performed in the framework of Einstein's 
theory of general relativity. For this purpose the broad collection of models for 
the equation of state, derived in part one, will be used as an input. Finally, in 
part three a comparison of the theoretically computed neutron star properties 
(e.g., masses, radii, moments of inertia, red and blueshifts, limiting rotational 
periods, cooling behavior) with the body of observed data will be performed. 
From this conclusions about the interior structure of neutron stars - and thus 
the behavior of matter at extreme densities - will be drawn. 

1 Introduction 

Neutron stars contain matter in one of the densest forms found in the universe. Matter 
in their cores possesses densities ranging from a few times po to an order of magnitude 
higher. Here po = 0.15 nucleons/fm 3 denotes the density of normal nuclear matter, 
which corresponds to a mass density of 2.5 x 10 14 g/cm 3 . The number of baryons 
forming a neutron star is of the order of A m 10 57 . The understanding of matter 
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under such extreme conditions of density is one of the central but also most complex 
problems of physics. 

Neutron stars are associated with two classes of astrophysical objects: Pulsars [[!]], 
which are generally accepted to be rotating neutron stars (the fastest so far observed 
ones have rotational periods of P = 1.6 ms, which corresponds to about 620 rotations 
per second), and compact X-ray sources (e.g., Her X-l and Vela X-l), certain of which 
are neutron stars in close binary orbits with an ordinary star. The first millisecond 
pulsar was discovered in 1982 ||, and in the next seven years about one a year has 
been found. The situation has changed radically with the recent discovery of an 
anomalously large population of millisecond pulsars in globular clusters 0, where 
the density of stars is roughly 1000 times that in the field of the galaxy and which are 
therefore very favorable environments for the formation of rapidly rotating pulsars 
that have been spun up by means of mass accretion from a binary companion. At 
present about 700 pulsars are known, and the discovery rate of new ones is rather 
high. 

As just outlined, neutron stars are objects of highly compressed matter so that 
the geometry of space-time is changed considerably from flat space. Thus for the 
construction of realistic models of rapidly rotating compact stars one has to resort to 
Einstein's theory of general relativity ||, |5|, || |7|, || || . The equation of state of the 
stellar matter, i.e., pressure as a function of energy density, is the basic input quantity 
whose knowledge over a broad range of densities (ranging from the density of iron 
at the star's surface up to ~ 15 times the density of normal nuclear matter reached 
in the cores of massive stars) is necessary in order to solve Einstein's equations. 
Unfortunately the physical behavior of matter under such extreme densities as in 
the cores of massive stars is rather uncertain and the associated equation of state is 
only poorly known. The models derived for it differ considerably with respect to the 
dependence of pressure on density, which has its origin in various sources. To mention 
several are: (1) the many-body technique used to determine the equation of state, (2) 
the model for the nucleon-nucleon interaction, (3) description of electrically charge 
neutral neutron star matter in terms of either (a) only neutrons, (b) neutrons and 
protons in (5 equilibrium with electrons and muons, or (c) nucleons, hyperons and 
more massive baryon states in f3 equilibrium with leptons, (4) inclusion of pion/kaon 
condensation, and (5) treatment of the transition of confined hadronic matter into 
quark matter. It is the purpose of this work to outline the present status of dense 
matter calculations, explore the compatibility of the properties of non-rotating and 
rotating compact star models, which are constructed for a collection of equations of 
state which accounts for items (l)-(5) from above, with observed data, and investigate 
the properties of hypothetical strange (quark matter) stars [[H| [TT], |l^]. The latter 
would constitute the true nature of neutron stars if strange quark matter should be 
more stable than ordinary nucleonic matter. 



2 Equation of state of dense neutron star matter 
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2.1 Theoretical framework 



2.1.1 Non-relativistic approach 

For non-relativistic models, the starting point is a phenomenological nucleon-nucleon 
interaction. In the case of the equations of state reported here, different two-nucleon 
potentials (denoted V^) which fit nucleon-nucleon scattering data and deuteron prop- 
erties have been employed. Most of these two-nucleon potentials are supplemented 
with three-nucleon interactions (denoted Vijk). The hamiltonian is of the form 

;2\ 
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(1) 



The many-body method adopted to solve the Schroedinger equation is based on the 
variational approach |TJ|, [TJ], [15] where a variational trial function | \P„ > is constructed 
from a symmetrized product of two-body correlation operators (Fy) acting on an 
unperturbed ground-state, i.e., 



s n f * 

i<j 



1$ > 



where |$ > denotes the antisymmetrized Fermi-gas wave function, 

|$ > = Af[ expiipj-Xj) . 



(2) 
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The correlation operator contains variational parameters which are varied to minimize 

(4) 



the energy per baryon for a given density p (see Refs. |13|, |Tj, |15, [16| for details): 
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As indicated, E v constitutes an upper bound to the ground-state energy E . The 
energy density e(p) and pressure P(p) are obtained from Eq. (f|) by 



e(p) = p \E v (p) + m) , 



P(p) 



9 d j . 
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(5) 



which leads to the equation of state in the form P(e) used for the star structure 
calculations here. 



2.1.2 Relativistic approach 

As discussed elsewhere H [L7|], the lagrangian governing the dynamics of many-baryon 
(-lepton) neutron-star matter has the following structure: 
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Table 1: Masses, electric charges (qs) and quantum numbers (spin (Jb), isospin (Is), 
strangeness (Sb), hypercharge (Y B ), third component of isospin (ias)) of the baryons 
included in the determination of the equation of state of neutron star matter. 



Baryon (B) 

■J \ / 


m B [Me VI 

1J L J 


Jb 


Ib 


S b 


Yb 


I3B 


q B 


fl 


939.6 


1/2 


1/2 

Lj Z 


o 


1 


-1/2 

±/ z 





D 


938.3 


1/2 

Lj Z 


1/2 

X/ Z 
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1/2 

Lj Z 


1 




1189 


1/2 


1 


-1 





1 


1 


s° 


1193 


1/2 


1 


-1 











Yr 


1197 


1/2 


1 


-1 





-1 


-1 


A 


1116 


1/2 





-1 











^0 


1315 


1/2 


1/2 


-2 


-1 


1/2 
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1/2 


-2 


-1 


-1/2 


-1 
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1232 


3/2 


3/2 





1 


3/2 
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A+ 


1232 


3/2 


3/2 





1 


1/2 


1 


A 


1232 


3/2 


3/2 





1 


-1/2 





A- 


1232 


3/2 


3/2 





1 


-3/2 


-1 



The subscript B runs over all baryon species that become populated in dense star 
matter, which are listed in Table [1|. The nuclear forces are mediated by that collec- 
tion of scalar, vector, and isovector mesons (M) that is used for the construction of 
relativistic one-boson-exchange potentials [|T^, |T^| . The equations of motion resulting 
from Eq. (^) for the various baryon and meson field operators (given in [|| [lTj ) are to 
be solved subject to the conditions of electric charge neutrality, pf ot = PBary+Piip = 0) 
which leads to 



j2q B (2j fl + i; 

B 



PF,B 

Qn 2 



E 



Pf,x 
3tt 2 







and (3 (chemical) equilibrium |20"1 , 



H B = (i n - q B H e 



(7) 



(8) 



where [i n and fi e denote the chemical potentials of neutrons and electrons, whose 
knowledge is sufficient to find the chemical potential of any other baryon of electric 
charge qs present in the system. To solve the equations of chemically equilibrated 
many baryon matter we apply the Greens function method, which is outlined next. 

The starting point is the Martin- Schwinger hierarchy of coupled Greens functions 
|2"2"| , |2"3"|| . In the lowest order, the Martin-Schwinger hierarchy can be truncated by 
factorizing the four-point Greens function ^(l) 2; 1'2') [unprimed (primed) arguments 
refer to ingoing (outgoing) particles] into a product of two-point Greens functions 
g [= gi(l'; 1')]. This leads to the well-known relativistic Hartree (i.e., mean- field) 
and Hartree-Fock approximations. The T-matrix approximation (also known as A 
or ladder approximation), which goes beyond this, truncates the Martin-Schwinger 
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hierarchy by factorizing the six-point Greens function g 3 (123; 1'2'3') into products of 
four- and two-point functions by which dynamical two-particle correlations in matter - 
which are connected with the two-body potential (denoted v ) - are taken into account. 
The main problem which one encounters hereby is the calculation of the effective 
scattering matrix (effective two-particle potential) in matter, T, which satisfies 

T = v — v ex + J v AT . (9) 

We have restricted ourselves to the so-called A 00 approximation, for which the nucleon- 
nucleon propagator is given by the product of two free two-point Greens functions, 
i.e., A = A 00 = ig°g°, and the relativistic Brueckner-Hartree-Fock (RBHF) approach, 
where both intermediate baryons feel the nuclear background medium and scatter only 
into states outside the Fermi sea [ZHj. The basic input quantity in Eq. ([]) is the 

nucleon-nucleon interaction in free space as derived, for example, in the Bonn meson- 
exchange model |19| . We have adopted the latest version of this interaction together 
with Brockmann's potentials A-C to compute the T matrix in neutron matter up to 
several times nuclear matter density []7|, [24j 25], p6| . The important feature of such 



meson-exchange models is that the potential parameters are adjusted to the two- 
body nucleon-nucleon scattering data and the properties of the deuteron, whereby 
(in this sense) a parameter-free treatment of the many-body problem is achieved. 
The "Born" approximation of T sums the various meson potentials of the nucleon- 
nucleon interaction in free space, i.e., 

< 12 | v | 1'2' >= £ 5{ v Tf v A^ r& , (10) 

M=a,Lu,ir,p,ri,S,<f) 

and thus neglects dynamical nucleon-nucleon correlations. It is this approximation 
which leads to the relativistic Hartree and Hartree-Fock approximations |], 27, 28, pT 



M 



The symbol T M in Eq. ( |10"D stands for the various meson-nucleon vertices, and A 
denotes the free meson propagator of a meson of type M. 

The nucleon self-energy (effective one-particle potential) is obtained from the T 



matrix as 0, 2£ 



= i£ / [tr (t bb> g B ') - T BB ' g B ' } , (11) 

B> 

where B' = p, n, S 1 * 1 ' , A, S ' - , A ++,+ ' 0,_ sums all the charged baryon states whose 
threshold densities are reached in the cores of neutron star models. The baryon 
propagators in Eq. ( pTTj) are given as the solutions of Dyson's equation, 

g B = g° B + g 0B Z B ({g B '}) g B , (12) 

which terminates the set of equations that are to be solved self-consistently. The 
diagrammatic representation of the second term in Dyson's equation, g 0B 'E B g B , is 
given in Fig. [I]. There, single lines denote the free propagator, g 0B , and double lines 
refer to the self-consistent propagator in matter, g B . This term corrects the free 
propagators (first term) for medium effects arising from the Fermi sea of filled baryon 
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Figure 1: Graphical representation of the self-energy contributions (second term in 
Eq. fll2f) ) arising from both the Fermi sea, i.e., nuclear matter consisting of filled 
baryon states of energies 0^(0) < u B (p) < fi B (upper graph), as well as Dirac sea 
(states — oo < uj B {p) < co> B (0), lower graph). The former (latter) lead to so-called 
medium (vacuum) polarization contributions BO, 31]. 



states, i.e., the nuclear- matter medium (upper shaded area in Fig. [I]), within which 
the baryons move.]] Obviously these vanish for baryons propagating in free space, for 
which Y, B = 0. Note that the determination of g B from Dyson's equation leads to a 
self-consistent treatment of the coupled matter equations (B|)-(0). 

The equation of state follows from the stress-energy density tensor, T^ u , of the 
system as 

E{p) = < Too > I P — m , where (13) 

= J^M^fWj^-grM. (14) 

The sum in the latter equation sums the contributions coming from the baryons and 
leptons (A = e, /i). The quantity C denotes the lagrangian given in Eq. (||) (see Ref. 
|2T| for details). The pressure is obtained from E(p) via Eq. (H). 



1 Contributions coming from the lower shaded area account for vacuum polarization corrections. 
The so-called no-sea approximation, which has been applied for the determination of most of the 



equations of state presented in Sect. 2.2, neglects such corrections. A critical discussion of the influ- 
ence of vacuum renormalization on the equation of state of high-density matter has been performed 
in Ref. Q . It was found that these have negligible influence on the equation of state up to densities 
of at least ten times normal nuclear matter density, provided the coupling constants are tightly con- 
strained by the saturation properties of nuclear matter. In the present paper, vacuum polarization 
contributions are contained in equations of state denoted G300 and GJ 00 (see Table 0). 
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Table 2: Nuclear equations of state applied for the construction of models of general 
relativistic (rotating) neutron star models. 



Label 



EOS 



Description (see text) Reference 



Relativistic field theoretical equations of state 



w 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 



G300 
HV 

pDCM2 
^180 
pDCM2 
^265 

^300 

^200 

^■Bonn "I" 

pDCMl 

^225 

pDCMl 

^B180 

HFV 



HV 



A* B HF + HFV 



H, K=300 
R,K=285 

Q,K=2Q5,B 1 / 4 = 180 

H, 7^=265 

H, 7T, K=300 

H,vr,K=200 

H,AT=186 

H,if=225 

Q,K=225, B l ' A = 180 

H,A,fT=376 

H,A,K=264 



17. 21 



13 
1 

i 

I 

w 

13 



34 



34 



21 



n m, mi 



Non-relativistic potential model equations of state 

12 BJ(I) H7A 

13 WFF(UVi4+TNI) NP, fsT=261 

14 FP(V M +TNI) N, 7^=240 

15 WFF(UVi4+UVII) NP, K=202 

16 WFF(AVi4+UVII) NP, K=209 

17 MS94 N,iT=234 



m, 0, n, n 



2.2 Models for the nuclear equation of state 

A representative collection of nuclear equations of state that are determined in the 
framework of non-relativistic Schroedinger theory and relativistic nuclear field theory 
is listed in Table |2|. A few of them are graphically shown in Fig. [| where the pressure 
is plotted as a function of energy density (in units of the density of normal nuclear 
matter, eo = 140 MeV/fm 3 ). This collection of equations of state has been applied 
for the construction of models of general relativistic (rotating) compact star models, 
which will be presented in Sect. f|. The specific properties of these equations of 
state are described in Table |], where the following abbreviations are used: N = pure 
neutron; NP = n, p, leptons; 7r = pion condensation; H = composed of n, p, hyperons 
(E 1 * 1 ' , A, S°'~), and leptons; A = Ai 232 -resonance; Q = quark hybrid composition, 
i.e., n, p, hyperons in equilibrium with u, d, s-quarks, leptons; K = incompressibility 
(in MeV); B 1 ^ = bag constant (in MeV). Not all equations of state of our collection 
account for neutron matter in full (3 equilibrium (i.e., entries 13-17). These models 
treat neutron star matter as being composed of only neutrons, or neutrons and protons 
in equilibrium with leptons, which is however not the ground-state of neutron star 
matter predicted by theory [T7|, |37|, [43| . As an example of such an equation of state, we 



exhibit the FP(V 14 + TNI) model in Fig. 0. The relativistic equations of state account 
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for all baryon states that become populated in dense star models constructed from 
them. As representative examples for the relativistic equations of state, we show the 
HV, HFV, G 30 o, and Ggfg 1 ^ 1 models in Fig. |2|. A special feature of the latter equation 
of state is that it also (as Ggfg 1 ^ 2 , which is not shown in Fig. §) accounts for the possible 
transition of baryon matter to quark matter. One clearly sees in Fig. |2| the softening 
of the equation of state, i.e., reduction of pressure for a given density, at e > (2 — 3) e 
which is caused by the onset of baryon population and/or the transition of baryon 
matter to quark matter. The stiffer behavior of HFV in comparison with HV at high 
densities has its origin in the exchange (Fock) contribution that is contained in the 
former equation of state. An inherent feature of the relativistic equations of state is 



that they do not violate causality, i.e., the velocity of sound, given by v s = c ^dP/de, 
is smaller than the velocity of light (c) at all densities, which is not the case for the 
non-relativistic models for the equation of state (cf. Table |3|). Among the latter only 
the WFF(UV 14 + TNI) equation of state does not violate causality up to densities 
relevant for the construction of models of neutron stars from it. 

The nuclear matter properties at saturation density related to our collection of 
equations of state are summarized in Table []. The listed quantities are: binding 
energy of normal nuclear matter at saturation density, E/A\ compression modulus, 
K; effective nucleon mass, M* (= m*/m, where m denotes the nucleon mass); asym- 
metry energy, a sy . With the exception of Ag° nn + HV and Ab™ f + HFV, the coupling 
constants of the relativistic equations of state are determined such that these satu- 
rate infinite nuclear matter at densities in the range 0.15 to 0.16 fm~ 3 for a binding 
energy per nucleon of about —16 MeV. For the Ag° nn + HV and A™ HF + HFV equa- 
tions of state the saturation properties are determined by respectively the relativistic 
Bonn and Brockmann meson-exchange models for the nucleon-nucleon interaction 
whose parameters are determined by the free nucleon-nucleon scattering problem and 
the properties of the deuteron (parameter-free treatment). The influence of dynami- 
cal two-particle correlations calculated from the scattering matrix leads for these two 
equations of state to a relatively soft behavior in the vicinity of the saturation density. 
This is indicated by the rather small compression moduli K related to these equations 
of state. All non-relativistic equations of state of our collection, which are determined 
in the framework of the variational method outlined in Sect. |2.1.1| , contain the impact 
of dynamical two-particle correlations in matter, too. The correlations are calculated 
for different hamiltonians. With the exception of BJ(I) and MS94, the calculations 
are performed for the Urbana and Argonne two- nucleon potentials Vu, UV14 flU 



and AV14 f45j, respectively, supplemented by different models for the three-nucleon 
interaction. These are the density-dependent three-nucleon interaction of Lagaris 
and Pandharipande, TNI fl46| , and the Urbana three-nucleon model, UVII [44j]. One 
sees that nuclear matter is underbound by ~ 4 MeV for two of these equations of 
state. The corresponding saturation densities are in the range of 0.17 to 0.19 fm~ 3 , 
thus nuclear matter saturates at somewhat too large densities for these equations of 
state. (The empirical saturation density is ~ 0.15 fm~ 3 |47]].) Equations of state la- 
beled 13 and 14 lead to binding energies and saturation densities that are in good 
agreement with the empirical values which has its origin in the density-dependent 
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Table 3: Nuclear matter properties of the equations of state used in this work. 



Label 


EOS 


E/A 


A) 


K 


M* 




e/eo 1 






[MeVl 

L J 


[f 

1 


m 3 1 

J 


[MeVl 


[MeVl 


[MeVl 




1 


P 

^300 


— 16.3 


o 


153 


300 


0.78 


32.5 




2 


HV 


-15.98 


n 


145 


285 


0.77 


36.8 




(J 


pDCM2 
^BISO 


—16 


n 


16 


265 


796 


32 5 




4 


pDCM2 
^265 


—16.0 


o 


16 


265 


0.796 


32.5 




5 


^300 


— 16.3 


o 


153 


300 


0.78 


32.5 




6 


^200 


-15.95 


o 


145 


200 


0.8 


36.8 




7 


A§? + HV 
il Bonn 


-11.9 


o 


134 


186 


0.79 






8 


pDCMl 
^225 


-16.0 


o 


16 


225 


0.796 


32.5 




9 


pDCMl 
^"BISO 


-16.0 





16 


225 


0.796 


32.5 




10 


HFV 


-15.54 





159 


376 


0.62 


30 




11 


A* B HF + HFV 


-14.81 





170 


264 


0.66 


32 




12 


BJ(I) 














23.1 


13 


WFF(UVi4+TNI) 


-16.6 





157 


261 


0.65 


30.8 


14 


14 


FP(Vi4+TNI) 


-16.00 





159 


240 


0.64 




5.6 


15 


WFF(UVi4+UVII) 


-11.5 





175 


202 


0.79 


29.3 


6.5 


16 


WFF(AVi4+UVII) 


-12.4 





194 


209 


0.66 


27.6 


7.2 


17 


MS94 


-16.04 





161 


234 




32.0 


13 



t Energy density in units of normal nuclear matter density beyond 
which the velocity of sound in neutron matter becomes larger (su- 
perluminal) than the velocity of light. The symbol "— " indicates that 
causality is not violated. 



three-nucleon interaction TNI.0 Interesting is the shift of the saturation density ob- 
tained for Ag° nn + HV and Ag BHF + HFV - relative to non-relativistic treatments that 
account for dynamical correlations too - toward smaller densities which is caused by 



relativity [|7], gg, |5§ |§ ]57| 



The four equations of state Gg c 5 M1 , Gg™ 2 , Gg™ 1 , and Gg^ 2 , which are based 
on the relativistic lagrangian of Zimanyi and Moszkowski, have only recently been 



determined |34], [58]. The transition of confined hadronic matter into quark matter 
is taken into account in equations of state Ggpg 1 ^ 1 (Fig. H) and Ggfg 1 ^ 2 . Here a bag 
constant of B 1 ^ = 180 MeV has been used for the determination of the transition 
of baryon matter into quark matter, which places the energy per baryon of strange 
matter at 1100 MeV, well above the energy per nucleon in 56 Fe (~ 930 MeV). Most 

2 It is well known that two-particle correlations alone fail in reproducing the empirical values of 
binding energy and saturation density. In this case the saturation points calculated from the standard 
Brueckner-Hartree-Fock jl3|, [l9| and non-relativistic T matrix approximations (so[ ||^] for 
different nucleon-nucleon interactions fall in a narrow band, often called the Coester band. It appears 
likely that this band would contain the calculated saturation point for any realistic nucleon-nucleon 
interaction 11531 Ml. 
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interestingly, the transition to quark matter sets in already at a density e = 2.3 eo 
|]33| , [34]], which lowers the pressure relative to confined hadronic matter. The mixed 
phase of baryons and quarks ends, i.e. the pure quark phase begins, at e ~ 15 eo, 
which is larger than the central density encountered in the maximum-mass star model 
constructed from this equation of state. We stress that these density thresholds are 
rather different from those computed by other authors in earlier investigations. The 
reason for this lies in the realization that the transition between confined hadronic 
matter and quark matter takes place subject to the conservation of baryon and elec- 
tric charge. Correspondingly, there are two chemical potentials, and the transition of 
baryon matter to quark matter is to be determined in three-space spanned by pres- 
sure and the chemical potentials of the electrons and neutrons. The only existing 
investigation which accounts for this properly has been performed by Glendenning 
p3| , |34j. Further important differences between the determination of Ggfg 1 ^ 1 and ear- 
lier (and thus inconsistent) treatments concern the description of the dense interior of 
compact starsQ and the approximation of the mixed phase as two components which 
are separately charge neutral. 



If the dense core may be converted to quark matter p9L K5C , 61 , it must be strange quark matter 



since 3-fiavor quark matter has a lower energy per baryon than 2-flavor. And just as is the case 
for the hyperon content of neutron stars, strangeness is not conserved on macroscopic time scales. 



Many of the earlier discussions |>9|, |6fj, |61|, |62|, |63|, |6J, |65[ |66|, |67[ have treated the neutron star as 
pure in neutrons, and the quark phase as consisting of the equivalent number of u and d quarks. 
However neither is pure neutron matter the ground state of a star nor is a mixture of — 2n u \ In 
fact it is a highly excited state, and will quickly weak decay to an approximate equal mixture of u, 
d, s quarks. 
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Figure 3: Distribution of observed neutron star masses [70 . 



3 Observed neutron star properties 

The global neutron star properties such as masses, rotational frequencies, radii, mo- 
ments of inertia, redshifts, etc. are known to be sensitive to the adopted microscopic 
model for the nucleon-nucleon interaction or, in other words, to the nuclear equation 



of state Thus, by means of comparing the theoretically determined values for 
these quantities with observed ones one may hope to learn about the physical be- 
havior of matter at super-nuclear densities. In the following we briefly summarize 
important star properties. 



3.1 Masses 

The gravitational mass is of special importance since it can be inferred directly from 
observations of X-ray binaries and binary pulsars (e.g., the Hulse- Taylor radio pulsar 
PSR 1913+16 IH). Rappaport and Joss were the first who deduced neutron star 



masses for six X-ray binaries [|H]| . A reexamination of these masses became possible 
owing to the improved determinations of orbital parameters JF2]. Improved values 
are shown in Fig. |. The objects 3U1700-37 and Cyg X-l are two non-pulsating 
X-ray binaries. It is expected that at least one of these objects, likely Cyg X-l, is 
a black hole. Very remarkable is the extremely accurately determined mass of the 
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Hulse- Taylor binary pulsar, PSR 1913+16, which is given by 1.444 ± 0.003 [BS . 

In summary, Fig. § suggests that the most probable value of neutron star masses, 
as derived from observations of binary X-ray pulsars, is close to about 1.4 M Q , and 
that the masses of individual neutron stars are likely to be in the range 1.1 S M/M & IS 



1.8 P| 



3.2 Rotational frequencies of fast pulsars 

The rotational periods of fast pulsars provide conditions on the equation of state when 
combined with the mass constraint [JT"3|| . As already mentioned in Sect. [I], the fastest 



so far observed pulsars have rotational periods of 1.6 ms (i.e., 620 Hz). The successful 
model for the nuclear equation of state, therefore, must account for rotational neutron 
star periods of at least P = 1.6 ms as well as masses that lie in range listed in Sect. 

n. 



3.3 Radii 

Direct radius determinations for neutron stars do not exist. However, combinations of 
data of 10 well-observed X-ray bursters with special theoretical assumptions lead Van 
Paradijs ]7J| to the conclusion that the emitting surface has a radius of about 8.5 km. 



This value, as pointed out in [f7ql , may be underestimated by a factor of two. Fujimoto 
and Taam ]76| derived from the observational data of the X-ray burst source MXB 
1636-536, under rather uncertain theoretical assumptions, a neutron star mass and 
radius of 1.45 M & and 10.3 km. An error analysis lead them to predicting mass and 
radius ranges of 1.28 to 1.65 M & and 9.1 to 11.3 km, respectively. When comparing 
these values with computed neutron star data, however, one should be aware of the 
fact that burster are suspected, but not known, to be neutron stars. 

3.4 Moment of inertia 

Another global neutron star property is the moment of inertia, I. Early estimates of 
the energy- loss rate from pulsars J7^] spanned a wide range of /, i.e., 7 x 10 43 < I < 
7 x 10 44 g/cm 3 . From the luminosity of the Crab nebula (~ 2 — 4 x 10 38 erg/sec), 
several authors have found a lower bound on the moment of inertia of the pulsar given 
by / > 4 - 8 x 10 44 g cm 2 M, fig M. 



3.5 Redshift 



Finally we mention the neutron star redshift, z. Liang [^1|] has considered the neutron 
star redshift data base provided by measurements of 7-ray burst redshifted annihi- 
lation lines in the range 300 — 511 keV. These bursts have widely been interpreted 
as gravitationally redshifted 511 keV e ± pair annihilation lines from the surfaces of 
neutron stars. From this he showed that there is tentative evidence (if the interpre- 
tation is correct) to support a neutron star redshift range of 0.2 < z < 0.5, with the 
highest concentration in the narrower range 0.25 < z < 0.35. A particular role plays 
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the source of the 1979 March 5 7-ray burst source, which has been identified with 
SNR N49 by its position. From the interpretation of its emission, which has a peak 
at ~ 430 keV, as the 511 keV e ± annihilation line |82], |83| the resulting gravitational 
redshift has a value of z = 0.23 ± 0.05. 



4 Properties of neutron star models 
4.1 Non-rotating star models 

The structure of spherical neutron stars is determined by the Oppenheimer-Volkoff 
equations p|, RH| , 



dP _ e(r) m(r) [1 + P(r)/e(r)] [1 + Aixr z P(r) /m(r)] 
dr r 2 [1 — 2m(r)/r] 



(15) 



which describe a compact stellar configuration in hydrostatic equilibrium. (Here we 
use units for which the gravitational constant and velocity of light are G = c = 1. 
Hence M & = 1.5 km.) The boundary condition reads P(r = 0) = P c = P(e c ), 
where e c denotes the energy density at the star's center, which constitutes the free 
parameter that is to be specified when solving the above differential equation for 
a given equation of state. The latter determines P c and the energy density for all 
pressure values P < P c . The pressure is to be computed out to that radial distance 
where P(r) = 0, which determines the star's radius R, i.e., P(r = R) = 0. The mass 
contained in a sphere of radius r (< R), denoted by m(r), follows from e(r) as 

r 

m (r) = 4tt J dr' r' 2 e(r') . (16) 


The star's total (gravitational) mass is given by M = m(R). 

Figure |4] exhibits the gravitational mass of non-rotating neutron stars as a function 
of central energy density for a sample of equations of state of Table |2|. Each star 
sequence is shown up to densities that are slightly larger than those of the maximum- 
mass star (indicated by tick marks) of each sequence. Stars beyond the mass peak 
are unstable against radial oscillations and thus cannot exist stably (collapse to black 
holes) in nature. One sees that all equations of state are able to support non-rotating 
neutron star models of gravitational masses M > M(PSR 1913 + 16). On the other 
hand, rather massive stars of say M > 2 M Q can only be obtained for those models 
of the equation of state that exhibit a rather stiff behavior at super nuclear densities 
(cf. Fig. §). The largest maximum-mass value, M = 2.2 M Q , is obtained for HFV, 
which is caused by the exchange term contained in that model. Knowledge of the 
maximum-mass value is of great importance for two reasons. Firstly, quite a few 
neutron star masses are known (Sect. |^), and the largest of these imposes a lower 
bound on the maximum-mass of a theoretical model. The current lower bound is 
about 1.56 M Q [neutron star 4U 0900-40 (= Vela X-l)], which, as we have just seen, 
does not set a too stringent constraint on the nuclear equation of state. The situation 
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Figure 4: Non-rotating neutron star 
mass as a function of central density. 
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Figure 5: Non-rotating neutron star 
mass as a function of redshift. 



could easily change if an accurate future determination of the mass of neutron star 
4U 0900-40 (see Fig. |) should result in a value that is close to its present upper bound 
of 1.98 M Q . In this case most of the equations of state of our collection would be ruled 
out. The second reason is that the maximum mass can be useful in identifying black 
hole candidates Ppl E7[ 153]. For example, if the mass of a compact companion of an 



optical star is determined to exceed the maximum mass of a neutron star it must be a 
black hole. Since the maximum mass of stable neutron stars in our theory is 2.2 M Q , 
compact companions being more massive than that value are predicted to be black 
holes. 

The neutron star mass as a function of gravitational redshift, defined as 

1 , (17) 



1 - 2M/R 



is shown in Fig. |5] for the same sample of equations of state as in Fig. |j. One sees that 
the maximum-mass stars have redshifts in the range 0.4 < z 0.8, depending on the 
softness (stiffness) of the equation of state. Neutron stars of typically M ~ 1.5 M© 
(e.g., PSR 1913+16) are predicted to have redshifts in the considerably narrower 
range 0.2 < z < 0.32. The solid rectangle covers masses and redshifts in the ranges 
of 1.30 < M/M Q < 1.65 and 0.25 < z < 0.35, respectively. As outlined in Sects, 
and p.5|, the former range has been determined from observational data of X-ray 



burst source MXB 1636-536 [76], while the latter is based on the neutron star redshift 
data base provided by measurements of gamma-ray burst pair annihilation lines |^TJ . 
(Note the remarks in Sects. |3.3j and |3.5| concerning the interpretation of these data.) 



From the redshift value of SNR N49 (if correct) we predict a neutron mass star of 
1.1 5 M/M & 5 1.6, which is consistent with the observed mass range given in Sect. 
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Figure 6: Radius as a function of 
redshift for a sample of equations of 
state of Table |2|. 
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Figure 7: Moment of inertia as a 
function of mass for a sample of 
equations of state of Table 0. 



|3.1| . The relativistic equations of state set a narrower mass limit for SNR N49 given 
by 1.4 < M/M Q < 1.6. 

Figure || displays the radius as a function of gravitational redshift. The solid dots 
refer to the maximum- mass star of each sequence. Of course, stars at their termination 
points possess the largest redshifts, and these become the smaller the lighter the stars. 
Under the assumption that the annihilation line interpretation is correct (Sect. |3.5| ), 
SNR N49 is predicted to have a radius in the range of 10 to 14 km. The relativistic 
equations of state lead to a narrower radii range, i.e., 12.5 to 14 km. In general, 
the non-relativistic equations of state lead to smaller radii for a given redshift. The 
reason for this lies in the relatively soft (stiff) behavior at low (high) nuclear densities 
of the non-relativistic equations of state, which is less pronounced for the relativistic 
Hartree and Hartree-Fock equations of state. (The softness/stiffness of the relativistic 
Brueckner-Hartree-Fock equation of state, however, is qualitatively similar to the one 
obtained for the non-relativistic models.) Small radius values of star models are 
important in order to achieve rapid rotation. For that reason star models constructed 
for the non-relativistic equations of state possess limiting rotational periods that are 
smaller than those obtained for the relativistic equations of state, as will be discussed 



in Sect. |4.2| ). However, because of causality violation of the non-relativistic equations 
of state at high nuclear densities (Sect. p.2|) , this trend may be an artifact. 
In Fig. |7| we show the moment of inertia of neutron stars, given by |89 



tt/2 R(0) 

/(«) = 4, / m J * e— e Ut(nL). ^ • (18) 
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as a function of gravitational mass. In Sect. 2!1 we have pointed out that, in general, 
the inclusion of baryon population in neutron star matter as well as the possible 
transition of confined hadronic matter to quark matter causes a softening of the 
equation of state, which leads to somewhat smaller star masses and radii. From the 
functional dependence of I on radius and mass, which is of the form / oc R 2 M, one 
expects a relative decrease of the moment of inertia of star models constructed for 
such equation of state. Of course, the general relativistic expression for the moment 
of inertia, given in Eq. (|I8| ), is much more complicated. It accounts for the dragging 
effect of the local inertial frames (frequency dependence u(r, 9)) and the curvature of 
space-time |p9| . Nevertheless the qualitative dependence of / on mass and radius as 
expressed in the classical expression remains valid |68|. Estimates for the upper and 
lower bounds on the moment of inertia of the Crab pulsar derived from the pulsar's 
energy loss rate (labeled Rud72), and the lower bound on the moment of inertia 
derived from the luminosity of the Crab nebula (labeled Crab) |78|, [79], BO] are shown 



in Fig. [fj for the purpose of comparison. (The arrows refer only to the value of /oab 
and not to its mass, which is not known.) 



4.2 Rotating star models 
4.2.1 Minimal rotational periods 

Figures |8| and exhibit the limiting rotational periods of compact stars, which is set by 
the gravitational radiation reaction-driven instability |9y, [H], It originates from 
counter-rotating surface vibrational modes, which at sufficiently high rotational star 
frequencies are dragged forward. In this case, gravitational radiation which inevitably 
accompanies the aspherical transport of matter does not damp the modes, but rather 
drives them |53|, Q . Viscosity plays the important role of damping such gravitational- 
wave radiation-reaction instabilities at a sufficiently reduced rotational frequency such 
that the viscous damping rate and power in gravity waves are comparable ||95|| . The 
instability modes are taken to have the dependence exp[iu m (fl)t + imcp — t/r m {Q)], 
where uj m is the frequency of the surface mode which depends on the angular velocity 
Q of the star, (f> denotes the azimuthal angle, and r m is the time scale for the mode 
which determines its growth or damping. The rotation frequency Q at which it 
changes sign is the critical frequency for the particular mode, m (=2,3,4,...). It is 
conveniently expressed as the frequency, denoted by that solves |K| [y refers to 
the viscosity dependence, see below) 



fit 



uj m (0) 



m 



a, 



i u m 



(19) 



where 



I2m(m-1) M 



(20) 



V 2m + 1 

is the frequency of the vibrational mode in a non-rotating star. The time scales for 
gravitational radiation-reaction r g<m , and for viscous damping time [0, T v>m , are 
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given by 



2 (m-l)[(2m + l)!!] 2 / 2m + 1 \ ro f R\ m+1 
9 ' m "' 3 (m + l)(m + 2) \2m{m-l)j \M J ' 1 ' 

R 2 1 

r "' m ~ (2m + 1) (m - 1) i/ ' ( } 

respectively. The shear viscosity is denoted by v. It depends on the temperature, T, 
of the star [u(T) oc T~ 2 ]. It is small in very hot (T »s 10 10 K) and therefore young 
stars and larger in cold ones. A characteristic feature of equations (|l9|)-(p2|) is that 

merely depends on radius and mass (R and M) of the spherical star model, which 
eases solving these equations considerably. 

The functions a m and 7 m contain information about the pulsation of the rotating 
star models and are difficult to determine @, j98| . A reasonable first step is to replace 
them by their corresponding Maclaurin spheroid functions a m and •y m PO, We 



therefore take a m (ft m ) and j m (£l m ) as calculated in Refs. |99|, |100|| for the oscillations 
of rapidly rotating inhomogeneous Newtonian stellar models (polytropic index n—1), 
and Ref. |K| for uniform-density Maclaurin spheroids (i.e., n=0), respectively. Man- 
agan has shown that depends much more strongly on the equation of state and 
the mass of the neutron star model [through u> m (0) and r gjm , see Eqs. ( p0| and (|2TD) 



than on the polytropic index assumed in calculating a m [ [L01| . 

Figure || shows the critical rotational periods, at which emission of gravity waves 
sets in in hot (T = 10 10 K) pulsars, newly born in supernova explosions. Figure |9| 
is the analog of Fig. |8|, but for old and therefore cold compact stars of temperature 
T = 10 6 K, like neutron stars in binary systems that are being spun up (and thereby 
reheated) by mass accretion from a companion. One sees that the limiting rotational 
periods P T (= 2vr/^) are the smaller the more massive (and thus the smaller the 
radius) the star model (cf. Fig. |6|). A comparison between Figs. |§| and |] shows that 
the instability periods are shifted toward smaller values the colder the star, due to 
the larger viscosity in such objects. Consequently, the instability modes of compact 
stars in binary systems are excited at smaller rotational periods than is the case for 
hot and newly born pulsars in supernovae.0 The dependence of P T on the equation 
of state is shown too in these figures. One sees that the lower limits on P T are set by 
the non-relativistic equation of state labeled 16 due to the small radii values obtained 
for the star models constructed from it. The relativistic models for the equation of 
state generally lead to larger rotational periods due to the somewhat larger radii of 
the associated star models. 

The rectangles in Figs. |] and |9| denoted "observed" cover both the range of ob- 
served neutron star masses, 1.1 5 M/M & IS 1.8 as well as observed pulsar periods, 



4 Sawyer 102] has found that the bulk viscosity of neutron star matter goes as the sixth power of 
the temperature, as compared with a T~ 2 dependence for the shear viscosity which is treated here. 
This means that at T > 10 9 K the bulk viscosity would dominate over the shear viscosity and thus 
damp the gravitational-wave instability. In this case the instability periods would be shifted toward 

which will be discussed below. 



values that are relatively close to the Kepler period Pk [10C, 103 
The latter sets an absolute limit on stable rotation because of mass shedding. Bounds on Pk are 
given in Table ||. 
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Figure 8: Gravitational radiation- 
reaction instability period P T versus 
mass for newly born stars of temper- 
ature T = 10 10 K pi . 



Figure 9: Gravitational radiation- 
reaction instability period P T ver- 
sus mass for old stars of temperature 
T = 10 6 K pi. 



i.e., P > 1.6 ms. One sees that even the most rapidly rotating pulsars so far observed 
have rotational periods larger that those at which gravity-wave emission sets in and 
thus can be understood as rotating neutron or hybrid stars.[] The observation of pul- 
sars possessing masses in the observed range but rotational periods that are smaller 
than say ~ 1 ms (depending on temperature and thus on the pulsar's history) would 
be in clear contradiction to our equations of state. Consequently the possible future 
observation of such pulsars cannot be reconciled - in the framework of our collection 
of models for superdense neutron star matter - with the interpretation of such objects 
as rapidly rotating neutron or hybrid stars. This conclusion is strengthened by the 
construction of neutron star models that are rotating at their Kepler periods, Pk, 
at which mass shedding at the star's equator sets in. Therefore this period sets an 
absolute limit on rapid rotation, which cannot be overcome by any rapidly rotating 
star. On the basis of neutron star models constructed from the selection of equations 
of state studied here, the smallest Kepler periods are found to be in the range of 
0.7 < Pk ~ 1 m s, depending on the softness of the equation of state[] (cf. Table £|). 

5 Depending on their composition, compact star models are denoted as neutron, hybrid, or strange 
stars. Neutron stars consist of protons, neutrons and more massive baryons in (3 equilibrium with 
leptons; hybrid stars are compact stars which, in addition to baryons, also contain quarks in their 
dense cores; hypothetical strange stars consist of (3-flavor) strange quark matter which, by hy- 
pothesis, may form the absolute ground-state of strongly interacting matter (see Sect. |] for more 
details). 

6 An investigation of this period for neutron stars and other compact objects that is performed 
without taking recourse to any particular models of dense matter (but derives the limit only on 
the general principles that: Einstein's equations describe stellar structure, matter is microscopically 
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P K is given by @, |, 



106 



ft S £ , with « K = w + £ + ^£ . (23) 

Note that Pk can only be obtained by means of solving Eq. (|23|) self-consistently in 
combination with Einstein's equation, 

n KX - l -g KX Tl = 8nT K \e,P(e)) , (24) 

and the equation of energy-momentum conservation, T kX .\ = 0, which renders the 
problem extremely complicated. The quantities TZ^, g KX , and 1Z denote the Ricci 
tensor, metric tensor, and Ricci scalar (scalar curvature), respectively. The depen- 
dence of the energy-momentum tensor T kX on pressure and energy density, P and e 
respectively, is indicated in Eq. (|24f) . The quantities u, v, and if) in Eq. ( ^3|) denote 
the frame dragging frequency of local inertial frames, and time- and space-like metric 
functions, respectively. 



4.2.2 Bounds on properties of rapidly rotating pulsars 

We restrict ourselves to the properties of a rapidly rotating pulsar model having a 
mass of M ~ 1.45 M , as supported by the evolutionary history of supermassive 
stars |75|. The bounds on its properties, whose knowledge is of great importance for 



the interpretation of fast pulsar, are summarized in Table f|. The listed properties 
are: period at which the gravitational radiation-reaction instability sets in, P T (in ms) 
with star temperature listed in parentheses; Kepler period, Pk (in ms); central energy 
density, e c (in units of the density of normal nuclear matter); moment of inertia, / 
(in g cm 2 ); redshifts of photons emitted at the star's equator in backward (zb) and 
forward (zp) direction, defined by 

1 /2 

and the redshift of photons emitted from the star's pole, 

Zp (Q) = e ~ u{n) - 1 . (26) 

(V denotes the star's velocity at the equator j|, All other quantities are as in Sect. 
[4.2. 1| .) According to Table |J newly born pulsars observed in supernova explosions 
can only rotate stably at periods 5cT 1 ms. Half-millisecond periods, for example, are 
completely excluded for pulsars made of baryon matter. Therefore, the possible future 
discovery of a single sub-millisecond pulsar, rotating with a period of say ~ 0.5 ms, 



stable, and causality is not violated) has only recently been performed by Glendenning |105| |. He 
establishes a lower bound for the minimum Kepler period for a M = 1.442 Mq neutron star of 
Pk = 0.33 ms. Of course the equation of state that nature has chosen need not be the one that 
allows stars to rotate most rapidly. 
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Table 4: Lower and upper bounds on the properties of pulsars with M 
calculated for the broad collection of equations of state of Table ^| [|108|| . 



1.45 M Q , 



P T (10 6 K) P T (10 10 K) P K e c /eo log/ 



upper bound 


1.1 


1.5 


1 


5 


45.19 


1.05 


-0.18 


0.45 


lower bound 


0.8 


1.1 


0.7 


2 


44.95 


0.59 


-0.21 


0.23 



would give a strong hint that such an object is a rotating strange star, not a neutron 
star, and that 3-flavor strange quark matter may be the true ground-state of the 
strong interaction, as pointed out by Glendenning | 107 |. Old pulsar of T = 10 6 K 
(and mass M ~ 1.45 M & ) cannot be spun up to stable rotational periods smaller 
than « 0.8 ms. Again, the two fastest yet observed pulsars, rotating at 1.6 ms, are 
compatible with the periods in Table |], provided their masses are larger than 1 M & 
|104|| . For the purpose of comparison the Kepler period, below which mass shedding 
at the star's equator sets in, is listed too. It may play a role in cold (hot) pulsar 
whose rotation is stabilized by its large shear (bulk) viscosity value (see footnote 4 
on Sawyers calculation of the viscosity in dense nuclear matter). 



5 Strange quark matter stars 
5.1 The strange matter hypothesis 

The hypothesis that strange quark matter may be the absolute ground state of the 
strong interaction (not 56 Fe) has been raised by Bodmer ||109|| , Witten []nj, and 
Terazawa | 110|| . If the hypothesis is true, then a separate class of compact stars could 



exist, which are called strange stars. They form a distinct and disconnected branch 
of compact stars and are not part of the continuum of equilibrium configurations that 
include white dwarfs and neutron stars. In principle both strange and neutron stars 
could exist. However if strange stars exist, the galaxy is likely to be contaminated 
by strange quark nuggets which would convert all neutron stars to strange stars 
|107| , |111| , |112j| . This in turn means that the objects known to astronomers as pulsars 
are probably rotating strange matter stars, not neutron matter stars as is usually 
assumed. Unfortunately the bulk properties of models of neutron and strange stars 
of masses that are typical for neutron stars, 1.1 IS M/M Q IS 1.8, are relatively similar 
and therefore do not allow the distinction between the two possible pictures. The 
situation changes however as regards the possibility of fast rotation of strange stars. 
This has its origin in the completely different mass-radius relations of neutron and 
strange stars (see Fig. |1(J) fl5"9"l| . As a consequence of this the entire family of strange 
stars can rotate rapidly - i.e., considerably below one millisecond -, not just those 
near the limit of gravitational collapse to a black hole as is the case for neutron stars. 
Furthermore, as I shall discuss below, the cooling history of neutron stars and strange 
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Figure 10: Radius function of 
mass of a non-rotating strange star 
with crust . 



Fi gure 11: The ratio -/crust /Aotai &s 
a function of star mass. Rotational 
frequencies are shown as a fraction 



of the Kepler frequency, f2 K [03 



stars may be quite different. 

5.2 Hadronic crust on strange stars and pulsar glitches 

At the present time there appears to be only one crucial astrophysical test of the 
strange-quark-matter hypothesis, and that is whether strange quark stars can give 
rise to the observed phenomena of pulsar glitches. In the crust quake model an oblate 
solid nuclear crust in its present shape slowly comes out of equilibrium with the forces 
acting on it as the rotational period changes, and fractures when the built up stress 
exceeds the sheer strength of the crust material. The period and rate of change of 
period slowly heal to the trend preceding the glitch as the coupling between crust 
and core re-establish their co-rotation. The existence of glitches may have a decisive 
impact on the question of whether strange matter is the ground state of the strong 
interaction. 

The only existing investigation which deals with the calculation of the thickness, 
mass and moment of inertia of the nuclear solid crust that can exist on the surface of 
a rotating, general relativistic strange quark star has only recently been performed by 
Glendenning and Weber |S3. Their calculated mass-radius relationship for strange 



stars with a nuclear crust, whose maximum density is the neutron drip density, is 
shown in Fig. |l^. (Free neutrons in the star cannot exist. These would be dissolved 
into quark matter as they gravitate into the strange core. Therefore the maximum 
density of the crust is strictly limited by neutron drip. This density is about 4.3 x 
10 11 g/cm .) Since the crust is bound by the gravitational interaction (and not 
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by confinement, which is the case for the strange matter core), the relationship is 
qualitatively similar to the one for neutron and hybrid stars, as can be seen from 
Fig. 13. The radius being largest for the lightest and smallest for the heaviest stars 



(indicated by the soid dot in Fig. [TOD in the sequence. Just as for neutron stars the 
relationship is not necessarily monotonic at intermediate masses. The radius of the 
strange quark core, denoted -Rdrip, is shown by the dashed line. (A value for the bag 
constant of B 1 ^ = 160 MeV for which 3-flavor strange matter is absolutely stable has 
been chosen. This choice represents weakly bound strange matter with an energy per 
baryon ~ 920 MeV, and thus corresponds to strange quark matter being absolutely 
bound with respect to 56 Fe). The radius of the strange quark core is proportional to 
M 1 / 3 which is typical for self-bound objects. This proportionality is only modified 
near that stellar mass where gravity terminates the stable sequence. The sequence of 
stars has a minimum mass of ~ 0.015 M & (radius of ~ 400 km) or about 15 Jupiter 
masses, which is considerably smaller than that of neutron star sequences, about 



0.1 M ||113|| . The low-mass strange stars may be of considerable importance since 
they may be difficult to detect and therefore may effectively hide baryonic matter. 
Furthermore, of interest to the subject of cooling of strange stars is the crust thickness 
of strange stars [114]. It ranges from ~ 400 km for stars at the lower mass limit to 
~ 12 km for stars of mass ~ 0.02 M Q , and is a fraction of a kilometer for the star at 



the maximum mass [8S 



The moment of inertia of the hadronic crust, J crU st, that can be carried by a strange 
star as a function of star mass for a sample of rotational frequencies of Q = Qk, 
and is shown in Fig. 11 . Because of the relatively small crust mass of the maximum- 
mass models of each sequence, the ratio / C rust//totai is smallest for them (solid dots in 
Fig. [11]). The less massive the strange star the larger its radius (Fig. [n]) and therefore 



the larger both / C mst as well as /total- The dependence of / CT ust and /total on M is such 
that their ratio I C rust //total is a monotonically decreasing function of M. One sees 
that there is only a slight difference between / crust for Q = and Q = 

Of considerable relevance for the question of whether strange stars can exhibit 
glitches in rotation frequency, one sees that /crust / /total varies between 10~ 3 and ~ 10~ 5 
at the maximum mass. If the angular momentum of the pulsar is conserved in the 
quake then the relative frequency change and moment of inertia change are equal, 



and one arrives at [89] 



^ = 4^ > = f ~ ( 1(r5 " 1(r3 ) f> with °</< 1 - ( 27 ) 

16 Jo i * 

Here Iq denotes the moment of inertia of that part of the star whose frequency is 
changed in the quake. It might be that of the crust only, or some fraction, or all of 
the star. The factor / in Eq. ( p7|) represents the fraction of the crustal moment of 
inertia that is altered in the quake, i.e., / = |A/|// crust . Since the observed glitches 
have relative frequency changes AQ/Q = (10~ 9 — 10~ 6 ), a change in the crustal 
moment of inertia of / 5 0.1 would cause a giant glitch even in the least favorable 
case (for more details, see [jS9|). Moreover, we find that the observed range of the 



fractional change in the spin-down rate, fi, is consistent with the crust having the 
small moment of inertia calculated and the quake involving only a small fraction / of 
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that, just as in Eq. For this purpose we write 



891 



AQ ASl/Cl \AI\ AQ/Q J crust 

7T " a^T" ~ Mi/ti f — > (1 ° tol0)/ ' (28) 



where use of Eq. (\FJ\) has been made. Equation fl28| ) yields a small / value, i.e., 
/ < (10~ 4 to lO^ 1 ), in agreement with / < 10 _1 established just above. Here 
measured values of the ratio (AQ/Q)/ (AQ/Q) ~ 1CT 6 to 10~ 4 for the Crab and Vela 
pulsars, respectively, have been used. 

5.3 Thermal evolution of neutron and strange stars 

The left panel of Fig. [E| shows a numerical simulation of the thermal evolution 
of neutron stars. The neutrino emission rates are determined by the modified and 
direct Urea processes, and the presence of a pion or kaon condensate. The baryons 
are treated as superfluid particles. Hence the neutrino emissivities are suppressed 
by an exponential factor of exp(— A/kT), where A is the width of the superfluid 
gap (see Ref. 115 1 for details). Due to the dependence of the direct Urea process 



and the onset of meson condensation on star mass, stars that are too light for these 
processes to occur (i.e., M < 1 M & ) are restricted to standard cooling via modified 
Urea. Enhanced cooling via the other three processes results in a sudden drop of 
the star's surface temperature after about 10 to 10 3 years after birth, depending on 
the thickness of the ionic crust. As one sees, agreement with the observed data is 
achieved only if different masses for the underlying pulsars are assumed. The right 
panel of Fig. |T2| shows cooling simulations of strange quark stars. The curves differ 
with respect to assumptions made about a possible superfluid behavior of the quarks. 
Because of the higher neutrino emission rate in non-superfluid quark matter, such 
quark stars cool most rapidly (as long as cooling is core dominated). In this case one 
does not get agreement with most of the observed pulsar data. The only exception is 
pulsar PSR 1929+10. Superfluidity among the quarks reduces the neutrino emission 
rate, which delays cooling ||115|| . This moves the cooling curves into the region where 
most of the observed data lie. 

Subject to the inherent uncertainties in the behavior of strange quark matter as 
well as superdense nuclear matter, at present it appears much too premature to draw 
any definitive conclusions about the true nature of observed pulsars. Nevertheless, 
should a continued future analysis in fact confirm a considerably faster cooling of 
strange stars relative to neutron stars, this would provide a definitive signature (to- 
gether with rapid rotation) for the identification of a strange star. Specifically, the 
prompt drop in temperature at the very early stages of a pulsar, say within the first 10 



to 50 years after its formation, could offer a good signature of strange stars ||114|| . This 
feature, provided it withstands a more rigorous analysis of the microscopic properties 
of quark matter, could become particularly interesting if continued observation of SN 
1987A would reveal the temperature of the possibly existing pulsar at its center. 
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Figure 12: Left panel: Cooling of neutron stars with pion (solid curves) or kaon con- 
densates (dotted curve). Right panel: Cooling of M — 1.8 M Q strange stars with 
crust. The cooling curves of lighter strange stars, e.g. M > 1M©, differ only in- 
significantly from those shown here. Three different assumptions about a possible 
superfluid behavior of strange quark matter are made: no superfluidity (solid), su- 
perfluidity of all three flavors (dotted), and superfluidity of up and down flavors only 
(dashed). The vertical bars denote luminosities of observed pulsars. 



6 Summary 

This work begins with an investigation of the properties of superdense neutron star 
matter. Various models for the associated equation of state of such matter are intro- 
duced and their characteristic features are discussed in great detail. These equations 
of state are then applied to the construction of models of non-rotating as well as ro- 
tating compact stars, and an investigation of their cooling behavior. The theoretically 
computed neutron star properties are compared with the body of observed data. 

The indication of this work is that the gravitational radiation-reaction driven in- 
stability in neutron stars sets a lower limit on stable rotation for massive neutron stars 
of P rs 0.8 ms. Lighter ones having typical pulsar masses of 1.45 M & are predicted to 
have stable rotational periods P > 1 ms. This finding may have very important impli- 
cations for the nature of any pulsar that is found to have a shorter period, say around 
P pa 0.5 ms. Since our representative collection of nuclear equations of state does 
not allow for rotation at such small periods, the interpretation of such objects as 
rapidly rotating neutron stars would fail. Such objects, however, can be understood 
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as rapidly rotating strange stars. The plausible ground-state state in that event is the 
deconfined phase of (3-flavor) strange quark matter. From the QCD energy scale this 
is as likely a ground-state as the confined phase of atomic matter. At the present time 
there appears to be only one crucial astrophysical test of the strange quark matter 
hypothesis, and that is whether strange quark stars can give rise to the observed phe- 
nomena of pulsar glitches. We demonstrate that the nuclear solid crust that can exist 
on the surface of a strange star can have a moment of inertia sufficiently large that 
a fractional change can account for the magnitude of pulsar glitches. Furthermore 
low-mass strange stars can have enormously large nuclear crusts (up to several thou- 
sand kilometers ||116| , |1171| ) which enables such objects to be possible hiding places of 



baryonic matter. Finally, due to the uncertainties in the behavior of superdense nu- 
clear as well as strange matter, we conclude that no definitive conclusions about the 
true nature (strange or conventional) of observed pulsars can be drawn from cooling 
simulations yet. As of yet, they could be made of strange quark matter as well as of 
conventional nuclear matter. 
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